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SUMMARY 


TITLE: The Tunable XUV Imager (TXI) Sounding Rocket Payload 
P.I.: Dr. Leon Golub, Smithsonian Astrophysical Observatory 

The TXI was flown successfully on 21 June 2001 (36.199 US). All systems functioned as planned 
and image data were acquired find sent to the ground. Unfortunately, due to a parachute failure 
the payload was destroyed. In this report we summarize results from the flight and provide detailed 
inf ormation on the high resolution X-ray imaging detector which was developed as part of the 
program. 
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1 Scientific, Technical & Man- 
agement 

1.1 Scientific Objectives 

In this section we describe the program of scien- 
tific studies to be carried out with the Tunable 
XUV Imager (TXI) payload. The instrumenta- 
tion is described in the Sec. 1.2. 

NASA’s Living With A Stax (LWS) and Sun- 
Earth Connections programs are implementing 
a core program of solar studies relevant to the 
interests of the U.S. in general and of NASA in 
particular. Some of the goals of these programs 
will be addressed by a U.S.-Japan collaboration 
on the Solar-B mission and the partially con- 
temporaneous Solar Dynamics Observatory and 
the STEREO Mission. However, recent stud- 
ies (discussed below) indicate that observations 
which are not part of the current plan are need. 
High spatial resolution observations in multiple 
wavelengths covering a broad temperature range 
with good temperature discrimination are cru- 
cial to further progress in understanding coronal 
heating and dynamics. We therefore feel that it 
is worthwhile to pursue rapid, low-cost missions 
which can be done within the scope of the NASA 
Suborbital Program, and which take advantage 
of recent technological progress in order to con- 
struct a new generation of instruments which of- 
fer the possibility of advancing LWS/SEC goals 
in a timely manner, especially in areas not in- 
cluded in the present mission complement. 

This proposal describes instrumentation 
which has been built and flown in an initial low- 
resolution test flight and which fits these criteria 
in a high-quality, low-cost package. By rebuild- 


ing the already-tested and qualified payload and 
bringing the full-resolution version to comple- 
tion we will make available within a short time- 
frame a new type of instrument whose design 
takes into account the most recent findings in 
coronal physics, with the aim of advancing not 
only LWS and SEC science objectives, but also 
coronal physics in general. 

Our group has, for the past two decades, 
been engaged in a comprehensive program of 
solar studies encompassing observations, mod- 
elling and theory. We rely on a small core group 
of in-house staff, while taking full advantage of 
collaborative arrangements with outside experts 
as needed for succesful completion of the scien- 
tific studies. A crucial part of our program is 
the availibility to the theoretical and modelling 
efforts of basic new observational data. This 
has in large part been accomplished by construc- 
tion and flight of new instrumentation - such as 
the NIXT rocket and the TRACE and Solar-B 
XRT telescopes - which can provide the neces- 
sary input to the modelling and the necessary 
constraints on theory. 

Below we provide a brief overview of prob- 
lem areas in coronal physics which we feel to be 
most pressing, which most need to be addressed 
in order to advance our understanding of coro- 
nal strucutre, stability and dynamics, and for 
which new types of observational capabilities are 
needed. We also discuss how the new data which 
we will obtain from flight of the TXI rocket pay- 
load will help to solve the problems which are 
now facing us. 
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1.1.1 Coronal Heating, Loop Fine Struc- 
ture and the Multithermal Corona 

The region in the Sun’s atmosphere between 
the photosphere and lower corona is one of a 
very rapid transition in physical parameters: 
the plasma /? (the ratio of gas pressure to mag- 
netic pressure) falls off from roughly unity in the 
photosphere to about 0.01 in the corona, while 
the magnetic filling factor increases by the same 
ratio. The temperature increases by more than 
two orders of magnitude, and the density de- 
creases by the same amount. This all happens 
in a layer from 2000 to 5000 km (Aschwanden et 
al. 2002) across. Observations indicate that this 
boundary layer is very dynamic and variable, 
the most dynamic region of the Sun. Its mor- 
phology is unknown: it is probably not a true 
layer, but a zone of separate magnetic structures 
with different plasma temperatures. Below this 
region, in the interior of the Sun, the generation 
of the magnetic field takes place, and above it, 
in the solar corona, rapid and intense release of 
magnetic field energy leads to the emission of 
gamma rays, hard X-rays, and energetic parti- 
cles up to the energy of cosmic ray particles. 
All the energy that heats the overlying corona 
and fuels the massive eruptions and flares must 
be transported through this region of the solar 
atmosphere. 

Theoretical models have focussed in recent 
years on the various ways in which energy may 
be transported to the corona, and there dissi- 
pated, through the mediation of magnetic fields. 
Virtually without exception, these models have 
in common the feature that the actual dissipa- 
tion of energy transported to the corona occurs 


in spatially highly localized regions, which may 
however either be truly localized or distributed 
throughout a larger volume. There is also obser- 
vational support for the presence of unresolved 
fine structure in the corona, based on studies 
using density-sensitive line ratios (e.g., Dere et 
al. 1987; Linford & Wolfson 1988; R. Thomas 
1995, personal communication), which typically 
yield filling factors of 0.01-0.001 for the transi- 
tion region and corona. 



Figure 1: Coronal fine structure at 1 arcsec res- 
olution. Active region seen in Fe ix/x with the 
TRACE imager, 19 May 1998. 

The new data from the TRACE satellite, as 
shown in Fig. 1, are beginning to resolve this 
substructure. The localized heating regions im- 
plied by these observations are far smaller than 
the coronal structures which have until now 
been observed, leading to the ancillary problem 
of having to explain how the energy which is 
released can be transported throughout the vol- 
umes which need to be heated. 
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A clue to the process by which coronal loops 
are heated may be provided by the new TRACE 
high resolution XUV data. Fig. 2 shows the 
brightening of a structure which, in lower reso- 
lution images, would have been called a “loop.” 
The right-hand frame shows that the loop is 
actually composed of numerous fine threads of 
bright material, and the two frames at left show 
the way in which this “loop” brightened: a cen- 
tral bright filament forms, which is then accom- 
panied by other bright parallel filaments; these 
brighten all along their lengths at once, rather 
than as a wave of brightness down from the 
top or up from the bottom. The broadening 
of the structure as it brightens is actually due 
to the successive formation of numerous bright 
filaments parallel to the original one, so that 
problem of cross-field transport may be avoided. 

We believe that these data indicate coronal 
heating via currents, which would explain the 
brightening all along the length of the bundle of 
“threads” (it can no longer be called a “loop”). 
The formation of additional bright threads par- 
allel to the first one may indicate heating by re- 
turn currents, which are thought to be necessary 
in current heating models or in models that in- 
voke a quasi-separatrix layer. The problem now 
becomes one of explaining why certain bundles 
of field lines become activated and how the en- 
ergy deposition within a restricted class of field 
lines takes place. It appears to be TRACE’S 
high spatial resolution combined with its abil- 
ity to isolate a narrow temperature range that 
makes these new observations possible. This ca- 
pability is shared by the TXI, but will now be 
extended to ionization states not reached by the 


SOHO EIT instrument or by TRACE. 

A number of recent studies (Yoshida et al. 
1995; Takeda et al. 1995; Walker et al. 1993) 
have emphasized the importance of thermally 
differentiated imaging of the coronal plasma as 
part of the program to understand the heating 
and structuring of the solar corona. It has long 
been known that very different structures are 
seen in ground-based eclipse data in, e.g., the 
‘red line’ (Fe X) compared with the ‘green line’ 
(Fe XIV) (see e.g., Thkeda et al. 1994). In soft 
x-ray observations, this point has been summa- 
rized by noting that the observed coronal loop 
structures may be viewed as relatively isolated 
mini-atmospheres (Rosner et al 1978), each with 
its own temperature and density structure. Be- 
cause of the marked difference in observed coro- 
nal parameters within small spatial distances, 
it appears that the coronal heating mechanism 
can have markedly different effectiveness in por- 
tions of the atmosphere which are otherwise in- 
distinguishable. The result is that data taken 
in emission lines of the same element at only 
slightly different formation temperatures can 
show very different structures. The point is 
explicitly demonstrated by the TRACE obser- 
vations, such as the comparison below of an 
Fe IX/X view of the corona (left-hand panel), 
with Fe XII (right-hand panel); the difference 
in temperature is only logT ~ 0.2. 

It is effectively impossible to determine tem- 
peratures by passband ratios if the same struc- 
tures are not seen in the two passbands, and 
if we do not observe enough temperature inter- 
vals to see all of the structures. This problem 
is solved by the TXI in the temperature range 
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Figure 2: A sequence of images taken in the TRACE 173A passband, showing a loop brightening 
event with multiple fine threads of emission. Time between frames is «10 min., FOV ~1X2 
arcmin. This passband includes both Fe IX and Fe X, which may cause some blurring. The new 
TXI instrument images each line separately. 


log T ~ 5.8— 6.4 by our ability to image all of the 
ionization stages from Fe IX through Fe XIV. 
The plasma has, in effect, nowhere to hide. 



Figure 3: Detail of TRACE Fe ix/x (left) and 
Fe XII (right) views of an active region, showing 
strong differences in coronal structure with only 
a small difference in temperature. 


1.1.2 Plasma Flows in the Corona 

The question of steady flows in coronal loops 
has been widely studied starting with Meyer 
and Schmidt (1968) (see also Glencross (1980), 
Cargill & Priest (1980), Noci (1981) and Lan- 
dini & Monsignori Fossi (1981)). Considering 
mass-conserving flows from one side of a loop 
to the other, the fundamental conclusion is that 
flows will arise whenever the loop is subjected 
to asymmetric heating or a pressure gradient 
from one footpoint to the other. For small pres- 
sure differences the flows remain subsonic, and 
the structure of the loops is not substantially 
changed by the presence of a flow. As the forc- 
ing is increased, the flow speed increases until 
the velocity reaches the sound speed at a point 
(the critical point) in the tube. The downstream 
velocity is then supersonic. 

The observations of flows in coronal loops 
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have been restricted to lines in the chromo- 
sphere, where the flows are necessarily small 
(Craig and McClymont, 1986). At coronal 
heights and coronal temperatures flows of sev- 
eral hundred km/s are expected. High-cadence 
timelapse movies from the TRACE satellite 
show persistent outflows at the sound speed in 
long coronal loops emanating from the vicin- 
ity of sunspots, while comparison with the 
SOHO/CDS provides spectroscopic evidence of 
flows (Winebarger etal. 2001). 

With the TXI we will be able to form spar 
tially resolved velocity maps of the entire corona 
at wavelengths such as 180A (Fe Xl) and 211 A 
(Fe XIV), which are particularly well-separated 
from other strong lines. This is accomplished 
by positioning the variable narrow passband so 
that it is alternately leading and trailing the 
position of the emission line. A line shift due 
to line-of-sight velocity then causes a change in 
strength in the two images whose sign depends 
on the direction of the velocity and whose mag- 
nitude depends on the magnitude of the veloc- 
ity; the difference between the two images then 
provides a measurement of the velocities. In this 
manner the amplitude, spatial distribution and 
temporal evolution of flows in the corona will 
be directly mapped for the first time. For the 
bright regions, where we expect 1% relative pho- 
tometric accuracy between frames, we resolve 
flows as small as 50 km/sec, comparable to the 
velocities seen in the TRACE movies. 

As an added bonus, we note that the 
monochromator mirrors operate near the Brew- 
ster angle, so that only the cr-polarized com- 
ponent is reflected. The polarization of the 


XUV emission can therefore be measured sim- 
ply by rotating the payload 90° around the line- 
of-sight. This is easily done with the pointing 
system available at White Sands. The difference 
image in this instance directly measures the po- 
larization, if any is present. 


Summary 

We have built and flight-tested a new sounding 
rocket payload which can provide the next step 
in observational capability needed to carry out 
a timely and innovative series of coronal studies, 
as we have described above. Because the pay- 
load was destroyed by a parachute failure be- 
fore our scientific objectives could fully be met, 
we are proposing to rebuild and continue flying 
the payload. The Tunable XUV Imager (TXI) 
will obtain non-overlapping, spectrally resolved 
images of the corona in all of the emission lines 
formed at temperatures from <1 MK to >3 MK, 
including - for the first time - velocity and po- 
larization images of the EUV corona. 

As has traditionally been the case with 
sounding rocket payloads, we view this instru- 
mentation as having both intrinsic scientific 
merit and as a test bed for future satellite in- 
struments. The proposed investigation is a valu- 
able complement to the SOHO and TRACE sci- 
ence, and it will provide a major addition to 
the LWS and SEC programs. The flight of the 
high-resolution TXI will take place near Solar 
minimu m, with a scientific program emphasiz- 
ing quiet coronal studies and bright points. 
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Figure 4: Image of the EUV corona in the 
Ni XVII line formed at 3 MK, obtained in the 
21 June 2001 flight. 

1.2 Technical Approach 

The Tunable XUV Imager (TXI) will provide 
high resolution spectrographic observations of 
the solar atmosphere from the transition region 
to active region temperatures. The scientific 
merits of these observations have been discussed 
above; here we provide a description of the in- 
strumentation which is being built to carry out 
this scientific program. 

The problems associated with observing the 
multi-temperature corona imply that the instru- 
mentation which is used must be able to: a) 
cover a wide range of temperatures, and b) iso- 
late a narrow temperature interval within the 
broad range. Both of these requirements are 
met by the Tunable XUV Imager, which pro- 
vides non-overlapping high resolution spectro- 
heliogram images in individual coronal emission 


lines. This is accomplished by using soft x-ray 
multilayer coatings, which provide high reflec- 
tivity at large reflectance angles which until re- 
cently were not accessible in the spectral region 
shortward of 300 A. 

The TXI consists of a broadband multi- 
layer coated imaging optic, which is an f/20 
on-axis sphere specially coated for a broad 
wavelength response by Dr. David Windt of 
Columbia University. Within the broad pass- 
band, a narrow wavelength slice is selected by 
a pair of narrow band analyzer elements; these 
are identically coated flat mirrors arranged as a 
Cowan-Golovchenko (Cowan 1983) double crys- 
tal monochromator (Figs. 5,7). 

The broadband response of the imaging el- 
ement, and the tunable analyzer bandwidth 
within that passband, are shown in Fig. 6; note 
that the mirror pair gives «1.6 times narrower 
bandwidth than either element alone, thereby 
providing spectral resolution adequate to sepa- 
rate the strong Fe lines in the TXI passband. 
The geometry of the Cowan-type monochroma- 
tor has the highly desirable property that, with 
a single drive element the mirrors move auto- 
matically in such a way that the entrance and 
exit beams remain fixed, so that there is no im- 
age motion in the focal plane during wavelength 
selection. The physical and optical parameters 
of the TXI are listed in Table 1. 

In order to meet the scientific goals of the 
investigation, we have selected the wavelength 
range A171 — 21lA, which includes strong spec- 
tral lines ranging from O V (log T = 5.4) 
through Fe XXTV (log T = 7.3). A listing of 
the strongest lines in the TXI range is provided 
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CCD Camera Primary Mirror 

Figure 5: TXI instrument layout, showing entrance filter, double-crystal monochromator, broad- 
band focussing optic and CCD camera. 


Bandpasses for TXI (Strongest Lines Indicated) 



Wavelength 

Figure 6: TXI broad and narrow passbands, with the most prominent XUV emission lines indi- 
cated. 

in Table 2.2. The count rates listed are for active imaging will be achievable. In addition, spectral 
regions and include all mirror reflectivities, po- separation substantially higher than the FWHM 
larization effects and filter transmissions, thus of the analyzer elements will be achieved by dif- 
representing true focal plane intensities. The ferential measurements taken with a AA ~lA. 
backside thinned CCD sensor used in our HRX Given the strengths of the lines involved, the 
camera has a DQE of «80% in the 40-70 eV available exposure times and the analyzer ele- 
energy range. ment bandwidth of AA ~4A, we estimate that 

all of the lines listed in Table 2.1 can be iso- 
The most prominent lines in the TXI spectral lat(i<i at lhe 1% separation level. The He II line 
range are those of Fh IX through Fe XIV. They „ ymk ^ be „ rongly Iedllced by the lhrK . 
are well-separated in wavelength, so that ‘clean’ 
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Physical & Optical Parameters of the TXI 

A Range: 

171-210A @ 10% of peak 
171-230 @ 1% of peak 

Analyzer AA: 

4 A 

Collecting Area: 

75 cm 2 geometric 
2.5 cm 2 eff. area @ peak 

Focal Length: 

2.2 meters 

Detector Resolution: 

2K x 2K 

Pixel Size: 

0.8 arcsec 

Exposure Times: 

1-3 secs. /wavelength 
10 sec @ selected As 


reflection system; however, any residual inten- 
sity can be calibrated by moving the analyzer 
mirrors to the second order position at ~152A. 
We note also that the presence of multiplets, 
such as the three closely-spaced lines of Fe XI 
near 193A, will not cause image smearing since 
there is no spatial dispersion in this instrument. 

A schematic diagram of the TXI is shown in 
Fig. 5. The monochromator portion contains 
a pair of mirrors which will operate in paral- 
lel at angles near 45°. At smaller near-grazing 
angles, the angular size of the sun would cause 
wavelength mixing which dilutes the monochro- 
mator resolution. At near-normal incidence set- 
tings, large mirror motions are needed to ac- 
complish wavelength changes and, more sub- 
tlely, the rapidly varying 7r-polarization reflec- 
tivity creates a difficult calibration problem. 
Mirror settings near the Brewster angle sacri- 
fice half of the reflectivity on the first reflection 
(and none thereafter) but provide stable polar- 
ization, clean spectral selection and mirror mo- 
tion of only ±12°. The as-built monochromator 
is shown in Fig. 7. 



Figure 7: Optical bench portion of the TXI in- 
strument, including entrance aperture, cooling 
block, CCD camera and XUV monochromator. 

In addition to the pure spectral imaging, the 
TXI can also provide line-of-sight velocity mea- 
surements by alternately positioning the steep 
positive and negative gradients of the multilayer 
response function at the wavelength of a given 
spectral line. The Doppler shift then causes 
the line to move up or down the reflectivity 
curve and the difference between the two im- 
ages is a Doppler image; the sum is equivalent 
to the on-band image. For a relative photomet- 
ric accuracy of 1% in the difference images, we 
will be able to measure line-of-sight velocities 
of 50 km/sec, which compares favorably with 
the ~100 km/sec motions we are seeing in the 
TRA CE data. 

We note also that the operational wavelength 
range of the TXI can easily be changed by 
switching out the mirror elements. Thus, on 
subsequent flights we have the option of moving 
to other spectral regions (e.g. A250-400A) in or- 
der to satisfy an extended set. of scientific objec- 
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Strongest Emission Lines Within the TXI Passband 


Ion 

A (A) 

logT 

Countrate (cts/pix/s) 

e /pix/s 

Fe IX 

171.08 

5.9 

360 

8700 

0 V 

172.17 

5.4 

44 


0 VI 

172.94, 173.08 

5.5 

30 

- 

Fe X 

174.53, 177.24 

6.0 

390, 240 

9200 

FeXI 

180.42 

6.1 

950 

2.2 xlO 4 

Si XI/Fe XII 

186.88 

6.2 

770 

1.7xl0 4 

FeXI 

188.22 

6.1 

670 

1.5 xlO 4 

Fe XXIV 

192.03 

7.3 

variable 

variable 

Fe XII 

192.40, 193.52, 195.13 

6.2 

1300, 680, 630 

3xl0 4 

Fe XIII 

202.04, 203.82 

6.2 

310, 290 

lxlO 4 

Fe XIV 

211.32 

6.3 

180 

3500 


tives. For mechanical reasons the tuning range 
of the TXI monochromator section is limited to 
a factor of two between lowest and highest wave- 
lengths in the flight version; other factors, such 
as locations of desirable spectral lines, availabil- 
ity of light-blocking filter materials, and multi- 
layer coating performance, also must be taken 
into account. 

1.3 Observing Program. 

The scientific program for the TXI has three 
main observational objectives: 

• Obtain high spatial resolution spectroheli- 
ogram images with good signal-to-noise in 
the series of iron lines from Fe IX at A171 
through Fe XIV at A21lA. Eight wavelength 
settings are needed and exposure times of 
0.1, 1 and 3 seconds are taken at each. 

• Obtain deep exposures for some weaker 
lines in order to extend the temperature 


range covered; exposures of 10 second du- 
ration will be taken at O V and O VI. 

• Obtain exposures at intermediate wave- 
lengths for spectral separation and Doppler 
shift measurements; eight 3-second expo- 
sures for this objective. 

The total available observing time for our 
previous rocket flights has averaged 320 seconds, 
with a minimum of «300 and a maximum of 
«350. The CCD camera can be read out while 
the monochromator is moving to the next po- 
sition, since both activities require about the 
same amount of time (4 secs). The total observ- 
ing program outlined above will require 240 secs, 
to complete, which is consistent with the typical 
min imum time above altitude for White Sands 
launches, and falls well within the expected du- 
ration of solar fine pointing. If estimated sun- 
acquired times are longer, we will obtain deep 
exposures at selected wavelengths, in order to 
image the corona beyond 1.5R©. 


<% 
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2 Differential Emission Mea- 
sure (DEM) Determination 

Extensive effort has been put into developing 
methods for reconstruction of the amount of 
emitting material in the corona as a function of 
temperature (the DEM) from observations such 
as those provided by the TXI instrument. These 
results have been presented at several confer- 
ences and a paper has been submitted for pub- 
lication (Golub etal. 2004). 

2.1 Two methods for DEM recon- 
struction 

The differential emission measure (DEM) is a 
basic observable of the coronal plasma in a re- 
gion of interest dV. The DEM is a measure of 
the amount of emitting plasma: 

JnldV = I DEM(T) d T , (1) 

where n c is the electron number density. The 
intensity d* observed in channel i is 

di = jDEM(T)n(T)dT, (2) 

where r»(T) combines the known information, 
i. e., the temperature response function for chan- 
nel i with the plasma emissivity function. The 
DEM may be found, in principle, by inverting 
Eq. 2, but that is an ill-conditioned problem if 
the r{(T) overlap in temperature coverage (as 
they do for the XRT and AIA). We discuss two 
approaches to solving this inversion problem: an 
iterative search method, and inversion methods. 

2.2 Iterative methods 

We wish to find the “best” fitting DEM for 
a given set of observations in several spectral 


channels. Here we consider a set of images taken 
of an active region (AR). We assume that the 
structure of the AR loops is unchanged during 
the time it takes for the instrument to accumu- 
late all of the images. Before processing, the set 
of images is co- aligned and reduced to physical 
units (e.£., DN/s). In this section, we discuss 
how we estimate the DEM in a given pixel. Our 
procedure produces an iterative least-squares fit 
to the observations using a DEM represented 
by a spline with evenly spaced knots in logT c 
space. For the set of observations (all i) in 
pixel p, we wish to solve the simultaneous equa- 
tions of Eq. 2. With the forward modeling 
approach, we assume a starting DEM, apply 
Eq. 2, and compare the predicted observations 
for each filter with the real observations. The 
optimal DEM spline is found using IDL mpfit 
routines from Craig B. Markwardt which search 
towards a minimum x 2 statistic by varying the 
spline knots. These routines use the non-linear 
least-squares method based on the MINPACK-1 
routines (http : //www . net lib . org/minpack/). 
This method, including uncertainty analysis, is 
described in more detail in Golub etal 2004. 

2.3 Direct inversion methods 

Another class of solution methods is comprised 
of inversion algorithms. Eq. 2 can be written in 
matrix form as 

Rxt DEMt = d% , (3) 

representing a simultaneous set of linear equa- 
tions, where DEMt is a vector of length iVr, and 
di is a vector of length TV*. For Nt = the 
system is determined and there is a unique solu- 
tion (assuming no degeneracies in the response 
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functions). The inversion is easily solvable with 
standard linear algebra techniques (e.g., lower- 
upper decomposition). However, unless the se- 
lection of temperature bins (implicit in the re- 
sponse functions) happens to be closely consis- 
tent with the characteristics of the “true” DEM 
function, the solution DEM is often unphysical. 
For example, solution functions often show sig- 
nificant negative DEM values. Furthermore, the 
constraint Nt = N{ may imply poor resolution 
in temperature. 

For Nt > the system in Eq. 3 is underde- 
termined, and there is an infinite family of solu- 
tions. In other words, the a posteriori informa ir 
tion (t.e., the data) is insufficient to constrain 
the model to a unique solution. In this case, 
there is freedom to constrain the solution with 
a priori information (t.e., physical constraints, 
such as assumptions about positivity or smooth- 
ness). Regularization methods turn Eq. 3 into a 
well-posed minimization problem with the addi- 
tion of such a constraint, by trading the quality 
of fit to the data against the quality of fit to the 
assumed constraint. (The amount of trade-off 
is typically adjusted such that the \ 2 fit to the 
data is ~ N{.) Note that this approach requires 
an explicit choice of a priori assumptions. 

Without making any a priori assumptions, 
the method of singular value decomposition 
(SVD) can be used to describe the family of 
solutions to Eq. 3. If iZ*r is considered as a 
linear map from the DEM space of Nt dimen- 
sions to the data space of Ni dimensions, then 
the SVD algorithm provides the Ni DEM or- 
thonormal basis functions which solve Eq. 3. 
This reduces the inverse problem to a square 


system of equations with a determinable prin- 
cipal solution. SVD also provides (Nt - Ni) 
DEM orthonormal null-space functions, which 
add to the principal solution to generate the in- 
finite family of solutions to Eq. 3. (Note that 
there is no physical significance inherent to the 
principal SVD solution, although it happens to 
be the solution vector of shortest norm.) 

2.4 Comparing DEM reconstruction 
methods 

We compared the iterative and SVD (inversion) 
methods to evaluate their pros and cons for im- 
plementation of an automated algorithm. A va- 
riety of nine input DEMs were designed and 
folded through the AIA response functions to 
produce simulated observations. Then each set 
of observations was solved for a reconstructed 
DEM using three methods: (a) “Simple” inver- 
sion (where Nt — Ni)] (b) SVD principal inver- 
sion; sind (c) the iterative approach. The “Sim- 
ple” inversion method was included as a baseline 
for comparisons, even though it was expected 
to produce relatively poor reconstructions. We 
used six of the AIA’s EUV channels (Ni — 6), 
and solved for DEMs with Nt = 19 temperature 
bins spanning the range logT = 5. 6-7.4 . We 
have quantified three measures on the relative 
value of the three solution methods considered. 
These are recorded in Ihble 1. The first quan- 
tity is runtime. For each set of observations, we 
calculated the average rimtime of the method for 
100 runs on a Sim Ultra 10 workstation, with a 
440 MHz processor and 1.0 GB of RAM. As ex- 
pected, the inversion methods are substantially 
faster than the iterative approach: r(iterative) 
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DEM model => 

a 

b 

C 

d 

e 

f 

g 

h 

i 

Median 

T&ble 1.1: Runtimes (milliseconds) 

“Simple” inversion 

0.055 

0.054 

0.053 

0.053 

0.055 

0.056 

0.055 

0.054 

0.054 

0.054 

SVD 

1.3 

1.4 

1.4 

1.4 

1.4 

1.4 

1.3 

1.3 

1.3 

1.4 

Iteration 

1800 

590 

750 

750 

570 

510 

1100 

520 

440 

590 

Table 1.2: \ 2 between observations from model and estimated DEM 

“Simple” inversion 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

SVD 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Iteration 

0.10 

0.11 

0.003 

4.3 

19.6 

0.038 

3.2 

0.009 

67.4 

0.11 

'Ib.ble 1.3: Cross-correlation coeff. between model and estimated DEM 

“Simple” inversion 

0.75 

0.063 

0.098 

0.54 

0.048 

0.069 

0.22 

-0.096 

0.053 

0.069 

SVD 

0.89 

0.95 

0.92 

0.67 

0.72 

0.69 

0.54 

0.99 

0.95 

0.89 

Iteration 

1.00 

1.00 

1.00 

0.56 

0.42 

0.64 

0.48 

1.00 

0.99 

0.99 


Table 1: Quantitatively comparing DEM reconstruction methods for different DEM models. Mod- 
els {a,b,c} are Gaussian distributions for cool, warm, and hot plasma. Models {d,e,f} are “double 
isothermal” peaks for cool, warm, and hot plasma. Model g is a cool isothermal peak plus a hot 
isothermal peak. Models {h,i} are multithermal. The median values are simply the medians for 
this arbitrary set of nine models. 


r 


i' 


12 




~ 400 x r(SVD) ~ 10 4 x r( “Simple” inversion). 

The second tabulated quantity is the \ 2 fit 
to the observations. Note that inversion meth- 
ods reproduce the observations perfectly. This 
is a mathematical consequence when solving a 
well-determined (or underdetermined) system of 
equations. Iteration converges to the first min- 
imum it reaches from the starting guess, which 
may only be a (nonzero) local minimum in \ 2 - 

The third tabulated quantity is the cross- 
correlation coefficient between the input and so- 
lution DEMs. The iterative method is shown 
to most accurately reconstruct the DEMs, typi- 
cally. “Simple” inversion does not have the tem- 
perature resolution to accurately reproduce the 
input DEMs. The SVD method does a fair job, 
considering that the principal solution is only 
one of a family of solutions, some of which might 
be much better. (In fact, since we are solving 
over the same temperature bins for which the 
input DEMs are defined, the input DEM is defi- 
nitely a member of the SVD family of solutions. 
Of course, this would not be expected when an- 
alyzing observations of coronal plasma.) From 
these results, it is apparent that the iterative 
approach does well at reproducing our sample 
of model DEMs, giving positive solutions. Its 
ability to reproduce the observations (in terms 
of the x 2 statistic) varies greatly. The inversion 
methods are much faster, which is a prominent 
concern when designing a DEM data product 
for, the ALA, which produces images with 
16+ Mpix at a cadence of 10 seconds. The SVD 
principal solution reproduces the observations 
perfectly (before error analysis), amd does ad- 
most as good a job at reconstructing the DEM 


models as the iterative method does, although it 
cam return negative components. There is free- 
dom to impose some physical constraunts on in- 
version methods (such as positivity) while re- 
tauning computation speeds, by allowing the x 2 
statistic to go to N{. 


2.5 Sixteen Million DEMs Every Ten 
Seconds 

Cadculating DEMs at several points in an image 
is straightforward for nearly amy method, how- 
ever the true power of large field of view DEM 
observations comes when the DEMs are avail- 
able for every image set. Say, for example, there 
is a coronad loop visible in a few aidjacent cham- 
nels, amd we wish to test the hypothesis that the 
loop is isothermad. With ALA standaird cawlence 
we cam follow the evolution of this loop from 
appearance to disappearance with a ten second 
cadence. If the loop lasts for 100 minutes we 
will have an observational dataset that consists 
of 600 sets of 6 EUV images to form DEMs from. 
For this project we would limit our amalysis to 
the few hundred pixels surrounding the loop, say 
1/4 million pixels/fraune totad. It is apparent 
that if the visible loop is in fact a physicadly 
distinct structure in the corona, its evolution 
will be different from the surrounding plasma. 
There will be “features” within the loop’s DEM 
that will evolve together; the rest of the DEM 
in the loop cam be attributed to changes adong 
the line of sight, not associated with the loop in 
question. An isothermad loop will have a DEM 
feature that is within the temperature resolution 
of the telescopes used (~ 0.2 log T c ). The ability 
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to calculate DEMs at every point along a loop 
(and all adjacent points, too) will provide an un- 
precedented set of observational constraints on 
coronal loop models and corona heating theo- 
ries. As shown in Ihble 1.1, the run times for 
SVD are about 400 times faster than the iter- 
ative method. However, the iterative method 
provides a better fit to the DEMs. The iter- 
ative model depends on an initial guess of the 
DEM; if we are clever in our choice of the first 
guess maybe the runtime for the iterative model 
can be improved. In the cases described above a 
flat curve was used as the initial guess. Can we 
speed up the reconstruction for a time series of 
mega-pixel images via a ‘‘pre-conditioned” guess 
that takes advantage of the shared DEM in- 
formation contained within mega-pixel images? 
Four methods are compared: Brute Force - our 
standard method; for each image, solve for the 
DEM of each pixel independently using a flat 
first guess. Single Average - for each image, av- 
erage all pixels and calculate a “global” DEM; 
use this as a first guess for the DEM reconstruc- 
tion of each pixel. Time Dependent - for the first 
image, solve for the DEM of each pixel indepen- 
dently using a flat first guess; use the resulting 
DEMs as a first guess for calculations in second 
image; iterate to last image. Time Averaging - 
for first image, use a global DEM as a first guess 
for DEM reconstruction for each pixel; use the 
resulting DEMs as a first guess for calculations 
in second image; iterate to last image. 

3 Electronic Systems 

The TXI rocket Electronic system has four 
prime functions: 


1: Control the operation of a E2V 42-40 CCD 
detector. This would include the external shut- 
ter for exposure and readout sequencing. 

2: Control the operation of the TXI monochro- 
mator for wavelength selection. 

3: Control of the detector and H/K data. 

4: Control of the command instruction. 

The electronic design unitizes an on-board 
single board computer (SBC). It is partitioned 
for prime science operations (detector image 
data) and for the command, housekeeping, and 
instrument operations. The SBC is a Kontron 
845D, running on a passive PCI Bus back plane. 
The SBC has two serial ports for co mman d and 
housekeeping telemetry, an Ethernet port for 
ground diagnostics, and a Compact Flash disk 
to increase boot-up speed. The flight software 
will be based upon existing code written in C++ 
and will run on a Windows NT operating sys- 
tem. The SBC will control the image acquisition 
of the Moses 42-40 CCD detector. 

After appropriate processing and formatting 
by a dedicated PCI frame grabber board (Moses 
#30259), the SBC utilizing a packet protocol 
with error-correcting capabilities will transfer 
images to the rocket telemetry interface. The 
software will run a pre-programmed image se- 
quence but will be able to accept real-time 
commands to change the sequence based upon 
quick look data. The data transmission is sup- 
plied by a NASA 10 Mbit downlink system for 
real time viewing during the flight. The data 
are also stored on-board for post-flight retrieval 
should data transmission fail, or have excessive 
dropouts. A GSE system will have four main 
functions: displaying of real-time thumbnails of 
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the data, displaying of full images as teleme- 
try constraints allow, displaying of housekeep- 
ing telemetry, and a command uplink to control 
image sequencing and exposure levels. Both the 
flight and GSE systems will require only mini- 
mal modifications to their original TXI designs. 

The TXI system block diagram (Fig. 9) 
shows the TXI electronic system. NASA sup- 
plies the command telemetry, data telemetry, 
telemetry interface, flight batteries, and video 
telemetry for the Kalpha camera. NASA will 
provide and service the flight batteries. A latch- 
ing changeover relay will switch the experiment 
between external and internal power; this will 
be under NASA control. The experiment’s bat- 
tery power will be applied at two minutes before 
launch. External power will be used for all test- 
ing, except for flight system tests. 

3.0.1 The High Resolution X-ray (HRX) 
Imaging Detector 

For the past decade our group at the Center 
for Astrophysics has been engaged in a program 
to design, build and fly high quality electronic 
imaging devices for soft x-ray and XUV ap- 
plications. We have successfully flown a num- 
ber of such devices in the NEXT rocket pro- 
gram, including intensified CCD cameras which 
have provided real-time images of the sun at He 
n A304A and Fe XII A193A. Under a NASA 
Grant, we developed the High Resolution X- 
ray (HRX) imaging detector, a CCD camera 
designed for flight and sensitive to both visible 
light and XUV wavelengths. This camera was 
flown succesfully in the previous TXI flight. For 
the past year we have been working with the 


DLR in Berlin, Germany to develop a larger- 
format camera, which will be available during 
the current calendar year for flight in the rebuilt 
TXI. 

The new HRX camera uses a back- 
illuminated E2V 42-40 CCD (with no anti- 
reflection coating). This is a 2Kx2K device, 
with 13.5 [i pixels. The DQE at the TXI wave- 
lengths is « 85% (as measured at SAO), since 
we are working below the silicon L-edge. This is 
substantially better than the efficiency of either 
photographic emulsions or of phosphor coated 
CCDs. The readout noise level with moderate 
cooling (-60° C) is 20 e“ per pixel at a read- 
out rate of 1,00,000 pixels/sec, so that a high- 
quality image requires an average exposure of 
only «250 e _ per pixel. The large-format CCD 
camera is used to provide the main TXI images, 
and a video-format camera provides a real-time 
Ha image. 

The flight version of the HRX camera is 
shown in Fig. 10. 

3.1 Affiliated Scientists 

As part of our overall approach, which has 
historically included extensive collaborative ar- 
rangements with scientists outside of the Smith- 
sonian, we have initiated specific collaborations 
for the work described in this proposal. Travel 
support for visitors is in most cases provided by 
internal SAO funding. 

The following list gives the names of the indi- 
viduals who will be working with us and a brief 
description of the role each will play in this pro- 
gram. 

Dr. G. Peres, Osservatorio di Palermo. Will 
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participate in data analysis with emphasis 
on modelling of coronal structures using the 
FLASH code, for which they are co-developers; 
he will also provide participation by members of 
the Osservatorio staff and students from the U. 
of Palermo. 

Prof. R. Rosner, U. of Chicago. Will provide 
expertise in implementation of the FLASH loop 
modelling code at SAO, will participate in the 
data analysis and will involve students from U. 
of Chicago in the analysis. 

Dr. B. Schmieder, IAP Medoc, will coordi- 
nate observations simultaneous with the TXI 
launch, at the newly-upgraded Pic du Midi Ob- 
servatory; science programs will focus on coro- 
nal dynamics on a target-of-opportunity basis. 

Dr. N.O. Weiss, DAMTP (Cambridge, UK). 
Will provide theoretical analysis of coronal data 
in the area of magneto-convection, studies of 
sunspot structure and dynamo-related global 
analysis; emphasis will be on the flows in 
penumbral loops seen at 10 6 K. 

3.2 Impact of This Work 

The Tunable XUV Imager (TXI) is a new con- 
cept in coronal imaging, which opens up new 
possibilities for scientific studies. For the first 
time we will be able to obtain high spatial res- 
olution imaging in isolated individual emission 
lines, while also being able to observe a large 
number of successive ionization stages of the 
same element; thus all of the coronal plasma 
within the covered temperature range will be 
observed. In addition, the same telescope can 
be used to obtain velocity maps and polarizar 
tion maps at any desired wavelength. 


The TXI permits us to address, among oth- 
ers, the following questions: 

- What is the distribution of temperature and 
emission measure in the corona, and how is it 
structured? 

- Where are plasma flows present in the corona, 
and what is their magnitude? 

- Where is XUV emission from the hot coronal 
plasma polarized, and by how much? 

The answers to these questions are directly 
relevant to determining the coronal heating 
mechanism(s), and to establishing the reasons 
for instability in coronal structures. We antic- 
ipate that the TXI concept will lead to fundar 
mental advances in our understanding of coronal 
physics, especially when used with a high resolu- 
tion, large format detector, and eventually from 
a platform capable of carrying out observations 
over an extended time period. 

The TXI will provide support and help in 
planning observing programs for several NASA 
missions, such as TRACE, Solar Probe and 
STEREO, since it observes coronal emission 
lines which are not available to the imaging in- 
struments on those satellites, and it provides im- 
ages for spectroscopic instruments. The types 
of emission measure, temperature, velocity and 
polarization data obtained from TXI are of po- 
tential usefulness to RHESSI, especially if we 
launch the sounding rocket payload as part of a 
flare- wait campaign, so that the TXI is observ- 
ing during a flare. 
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Figure 8: The number of independent channels used in the DEM reconstrcution is crucial: reducing 
the number from 6 to 5 introduces are large degradation in the quality of the reconstruction. 



Figure 9: Schematic block diagram of the TXI electronics system. 
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Figure 10: The HRX camera utilizes an efficient compact design for flight applications. 
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